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Section I 

INTRODUCTION 


In Space Shuttle vehicle ground wind related analyses peak wind speeds 
are frequently more useful than the mean values found in the climatological 
records. Also, it is much easier to use peak speeds in the wind monitoring 
efforts that take place during the vehicle on-pad pre-launch activities. 

Extracting peak values from original ground wind records was found to 
be very expensive, even when the data were readily available. To avoid the 
expense and to circumvent the missing data a search was undertaken to find 
a procedure whereby satisfactory estimates of hourly peak speeds could be 
obtained from standard observed hourly wind measurements which are archived 
with the National Climatic Center. 

Wind gustiness characteristics at an anemometer site are largely determined 
by the roughness of the underlying terrain, the air mass stability, ami the 
height of the instrument. In order to estimate hourly peak winds, recent 
turbulence and stability theory has been applied to hourly wind observations 
u, producing such quantities as the friction velocity u AO » the Monin- 
Obukhov stability length L^, and the standard deviation of wind speed ^ . These 
quantities lead to an estimate of peak wind U* , and there is an opportunity 
to optimize factor F in the final expression, u^ * u + F o 

To test the effectiveness of the theory on Cape Kennedy wind data, 

2952 observations of u and u were selected from a set of 12 midseasonal 

P 

hurricane-free months. As a control procedure, a simple optimized gust 
factor G was found from the same data sample and applied in the relation, 

Up * Gu. Then identical statistical evaluations were performed on peak wind 
estimates from all three procedures. 
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Section II 
PROCEDURE 

Two closely related methods were tested in the estimation of peak 10 m 
winds at Kennedy Spacecraft Center (KSC). The results of these tests on a 12- 
month sample were compared to a set of control statistics for the same 
sample. 

The general procedure follows some suggestions made by Dr. George H. 
Fichtl, NASA, MSFC. In brief, the preliminary objective is to find a mixing 
length estimate - the Monin-Obukhov stability length - from consideration of 
the Pasquill stability class: 

The solar radiation (I) is computed from the sun’s hour angle (t), its 
declination angle (6), its zenith angle (6), and the latitude of KSC ($) by 
equation 1, 

I » [1.585 cos 6] [1 - 0.19 (sec 0) 1 ' 2 ] (1) 

where 

cos 6 ■ cos t cos 6 cos <p + sin 6 sin <p 

The corrected solar radiation (I c ) Is found from equation 2, 

I - 1(1 - 0.01 aC) (2) 

c 

where 


a - 0.3 

c ■ cloud cover in percent . 

This correction is insensitive to cloud type, the smallest possible trans- 
mission being 0.7. Since transmissions in dense low overcasts can be as 
small as 0.2, an attempt was made to take account of individual cloud type, 
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height, and amount. A search of the literature revealed nine reports in 
Russian which are not readily available and papers by Lumb (ref. 1) and 
Parker (ref. 2). Lumb's paper presents an empirical method of estimating 
total solar radiation, based upon data taken over the ocean at 52 1/2°N 
latitude, and Parker's paper does the same thing for an island station near 
the Equator. However, cloudiness at the former site is judged to be too 
unlike that at KSC to apply Lumb's formulae, and Parker has utilized certain 
data which are not contained in KSC surface reports. Therefore the 
parameter "a" has been held constant throughout a particular trial run. 

The Pasquill stability class evaluation takes into account the 10-meter 
wind speed, the corrected solar radiation, and the nighttime cloud coverage. 
Table 2-1 and Figure 2-1, both from Pasquill (ref. 3), reveal six categories 
ranging from "extremely unstable" to "moderately stable". Examples of the 
utilization of a Pasquill stability classification may be found in estimates 
of plume dispersion, as by Turner (ref. 4). 


Table 2-1. PASQUILL STABILITY CLASSES 


A - Extremely unstable conditions 

D - Neutral conditions* 


B - Moderately unstable conditions 

E - Slightly stable conditions 

C - Slightly unstable conditions 

F - Moderately stable conditions 





NIGHTTIME CONDITIONS 

SURFACE WIND 

DAYTIME INSOLATION 

THIN OVERCAST 
OR >4/8 
CLOUDINESS 1 ' 

v3/8 

CLOUDINESS 

SPEED, M/SEC 

STRONG 

MODERATE 

SLIGHT 

<2 

A 

A-B 

B 



2 

A-B 

B 

C 

E 

F 

4 

B 

B-C 

C 

D 

E 

6 

C 

C-D 


D 

D 

>6 

C 

D 

D 

D 

D 

4 Applicable to heavy overcast t day or night 



i'The degree of cloudiness is defined as that fraction of the sky above 1 

the local apparent horizon which is covered by clouds. 
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The Monin-Obukhov stability length is then obtained by computerizing 
Table 2-1 and Figure 2-1. The roughness length it assumed to be 0.2 m 
at KSC. The subsequent step marks the point of divergence for the two 
trial procedures and the control procedure, which will be designated A, B, and 
C, respectively* 


Method A * The friction velocity u^ Q Is computed by equation 3, following 


Panofsky (ref. 5) and Paulson (ref. 6), 


u *o * ku/[ln(z/* o ) - i|»(z/L o )J 


(3) 


k - 0.4 


u * hourly wind speed, an approximation to the mean wind speed over the 
60-minute period of observation of peak winds 

z/L q ■ stability parameter 


z ■ 10 a 


If z/L 


( < 0, x ■ (1 - 16z/L and ifr » In ~ 

1 o o 1-* 

- 2 tan 1 + f 

o 

• 0, i> m 0 
> 0, i l> m - 5z/L • 


Then the standard deviation of wind speed is estimated by equation 4, 


° u * 2 * 5 u *o ‘ 


(4) 


Finally, the peak wind speed for a 60-minute period centered upon the 
hourly observation is estimated by equation 5, 


u • u + F o 
P u 

F ■ factor to be optimized. 


(5) 
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Method B > The friction velocity is computed as in Method A, but the 
limitation that o u /u^ 0 * 2.5 is removed by introducing a quantity Ri, S) 
which was tested by Panofsky et al., (ref. 7). 


Ri ■ Richardson number 


S * stability parameter 


B - Ri S / 


m f-- * 

I o 




o 


<_ 0, Ri 

I > 0, Ri 


and S ■ 1 

o 



-1/4 




and S 


1 + 5 - 2 r 

o 


as may be seen in Figure 2-2 (from reference 7), the ratio o u /° AO takes on 
the following values as 2 varies: 


c /u._ - 2.75 when 100 B < - C.175 
|| 

aj u* Q - 2.40 - 1.26 (100 B) when - 0.175 < 100 B < 0.254 
0 u /i^ o ■ 2.38 when 100 B > 0.254 • 

Equation 5 is again used to find u^. 

Method C. The function <i>(z/L ) is set equal to zero except in the case 

“ O 

where daytime Insolation is strong and the hourly observed wind exceeds 6 m/sec. 
This implies that neutral conditions are assumed to exist in most of the cases. 


In the computation of the peak wind, the value of u is multiplied by a 

gust factor G, where G ■ 1 + 0.16 F 1 and F' is optimized. The expression for 

piedicting u is then u ' ■ G u. 

P P 

In all of the above methods t the quantity u is really a short-period 
estimate of the mean wind through the 6(^-mlnute period for which peak wind 
speed is to be predicted. 
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Section III 
DATA 


The three methods were tested on three years of midseasonal months' 
data for Cape Kennedy, namely, January, 1966-1968; April, 1965-1967; July, 
1965-1967; and October, 1965-1967. Eight reports per day at 01, 04, 07, 10, 

13, 16, 19, and 22 EST were processed. There were no hurricanes in the vicinity 
of KSC during these months. 

The recording anemometer providing the data for this study was a three- 
b laded Bendix-Frlez installed at a height of 10 m. This anemometer had a 
threshold value of about 2.5 kt, so each recorded zero value of hourly wind speed 
was questioned. All zero values were finally replaced by an estimated 
average subthreshold value of 1.1 kt. In addition to these replacements, two 
obvious mistakes in data copying were corrected by interpolated values. 

The distribution of the 2952 hourly observations of wind s^eed at Cape 
Kennedy is given in Figure 3-1. There is a skewed peak in the observations 
at 4 kt and no reported winds exceed 22 kt. A bias favoring even-numbered values 
over o * 1- numbered values is apparent over about half of the range, and it 
introduces an irregularity in the estimated peak winds which is permitted to 
remain. 
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Figure 3-1. FREQUENCY DISTRIBUTION OF HOURLY WIND SPEED AT CAPE KENNEDY FOR 
JANUARY, APRIL, JULY, AND OCTOBER IN THE PERIOD FROM 1 APRIL 1965 
TO 31 JANUARY 1968. ONLY HOURS 01, 04, 07, 10, 13, 16, 19, AND 22 
EST ARE INCLUDED 
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Section IV 
RESULTS 


Taking all 12 months' data at a time, optimum F and F' factors were 
found from calculations of minimum rms errors- The outcome is shown in 
Figure 4-1, where F for Method A is 2.95, F for Method B is 2.7, F' for Method 
C (not shown) is 3.8. The spread of minimum rms error is but 0.34 kt, with the 
order of rank being C, A, B. The failure of Method B to surpass the other less 
complex methods indicates that Figure 2-2 may not be as representative of 
KSC as of the other sites to which it pertains. 

The distributions of observed peak winds and estimates were placed on 
bar graphs (Figures 4-2, 4-3, 4-4). These figures all show a large surplus of 
estimates in the range 0.0 - 3.5 kt, a feature which may be attributed to 
anemometer inertia because this class contains all of the adjusted zero values 
mentioned in Section III. 

For higher peak winds, each method has a notable excess in frequency 
somewhere within the 3.5 - 11.5 kt range. The excess for Method A is the least, 
and it occurs at 5.5 - 7.5 kt. The excess for Method C is the greatest, and it 
occurs at 9.5 - 11.5 kt. The excess for Method B is at 3.5 - 5.5 kt. All 
methods have notable deficits in frequency within the 10 - 15 kt range with 
another excess just under 20 kt. 

The chi-square test is rendered inapplicable by the unevenness in the 
frequency distribution of observed peak wind speed, which is an empirical 
estimate of the true distribution of this quantity. However, as a tool for 
comparison of the three methods, chi-square values were computed. By this 
criterion, the order of rank (with chi-square values in parentheses) is A 
(288 kt), B (467 kt), C (610 kt). 

Further testing was conducted by use of the Kolmogorov-Smirnov (K-S) 
test. The null hypothesis H q is, "for identical observational periods, there 
is no significant difference between the distributions of estimated peak wind 


4-1 


ERROR (KT 


— METHOD A 

METHOD B 

12 MONTHS DATA 


z o « 0.2 M 
a = 0.7 



Figure 4-1 . OPTIMIZATION CURVES FOR THE F- FACTOR OF METHODS A AND B 
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FREQUENCY DISTRIBUTION OF REPORTS OF PEAK WIND SPEED AND ESTIMATES 
OF PEAK WIND SPEED BY METHOD A. THE PERIOD IS AS NOTED IN 
FIGURE 3-1 
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Figure 4-3. FREQUENCY DISTRIBUTION OF REPORTS OF PEAK WIND SPEED AND ESTIMATES 
OF PEAK WIND SPEED BY METHOD B. THE PERIOD IS AS NOTED IN FIGURE 
3-1 
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speeds". Cumulative percent frequency (CPF) curves for the K.-S test are given 

in Figure 4-5, and the results are listed in Table 4-1. According to tabulated 

critical values of the K-S statistic, H need not be rejected at the five 

o 

percent level nor the one percent level. The ranking of the three procedures 
by this criterion is A, C, B. 

Table 4-1. KOLMOGOROV-SMIRNOV TEST VALUES FOR METHODS A, B, AND C 


A 1 

_B 

_N JC 

_N 

.086 14 

.112 

14 .103 

14 

N = number of degrees of freedom 



An additional comparison is available in a case by case tally of the most 
accurate method. This was carried on for the entire sample and for the subset 
of peak winds greater than 20 kt. Table 4-2 gives the outcome of the two 
counts, and it discloses an apparent superiority in Method C, the control 
procedure. This is especially evident in the subset of greater wind speeds, 
where Method C is best in 62.8 percent of the cases. 

Table 4-2. RELATIVE ACCURACY OF ESTIMATES OF PEAK WIND SPEED BY 
METHODS A, B, AND C. I = INDETERMINATE, E.G., THE 
EDGE IN ACCURACY IS LESS THAN 0.1 KT. 





METHOD 


All 

cases, frequency 

_A_ 

261 

B 

622 

C 

1582 

I 

489 

u p 

> 20 kt, frequency 

39 

49 

157 

5 

All 

cases, frequency (percent) 

8.9 

21.0 

53.6 

16.5 

U P 

> 20 kt, frequency (percent) 

15.6 

19.6 

62.8 

2.0 


The correction factor, a, of equation 2 was also set equal to a greater 
value of 0.7 to test the sensitivity of the procedures. The overall result 
was a reduction of rms error of 0.05 kt and a small improvement in distribution 
of estimates. It is therefore concluded that 0.7 is a slightly more realistic 
correction factor than 0.3 at Cape Kennedy. When cloudiness reaches 100 
percent, the transmission of insolation to the ground is then 30 percent and 
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CUMULATIVE FREQUENCY CURVES FOR REPORTS OF PEAK WIND AND FOR 
ESTIMATES OF PEAK WIND BY METHODS A, B, AND C. THE DATA ARE 
NOTED IN FIGURE 3-1 




the depletion is 70 percent- These values are generally representative of low 
stratified overcasts. 

Individual estimates of peak wind are sensitive to the selection of the 
factor, a, in those infrequent cases when the Pasquill stability class (Table 2-1) 
is thereby changed. Crossing over a class boundary can vary the peak wind 
estimate by 10-15 percent. 

Figure 4-6 gives a further comparison of the three methods, using the 
absolute error in peak wind as a criterion. The ranking is C, A, B, with this 
order being set in the 0-1 kt error class and continued throughout. The graph 
also discloses that about 30 percent of all cases are under 1 kt error, 80 
percent are under 3 kt, and 95 percent are under 6 kt error. 

Figure 4-7 uses percentage error in peak wind as a criterion, and it 
points up the lack of clear superiority of any method through the entire range 
of wind speeds. For examnle, at 5 kt, which includes observed hourly wind 
values from 3 to 7 kt, 90 percent of the cases of Method B are within 33 per- 
cent in accuracy, but this superiority does not hold through all of the per- 
centile curves at 5 kt. Perhaps the most interesting feature of the three 
graphs is the dip in the control procedure. Method C, at 10 kt, which includes 
values from 8 to 12 kt, in the 90 percentile curve. Ninety percent of all 
cases in this category are under 20 percent error and this gives Method C its 
most notable advantage. However, the advantage is not retained at higher wind 
speeds. 

The roughness length, z q , was changed from 0.2 m to 0.1 m to test the 
sensitivity of all procedures to this parameter. The results indicated little 
sensitivity at any point, that is, outcomes were practically identical to the 
previous runs. 
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Figure 4-6. CUMULATIVE FREQUENCY CURVES FOR ABSOLUTE ERROR IN ESTIMATION OF 
PEAK WINDS BY METHODS A, B, AND C. THE DATA ARE AS NOTED IN 
FIGURE 3-1 








Section V 

CONCLUSIONS 


The advantages and disadvantages of the three evaluation schemes for 
peak wind are summarized in the following paragraphs. 

Method A: uj, ■ u + 2.5 F u* G , where u* Q ■ ku/(A-0(z/L o )) . 

The rms error is 3.01 kt over a 12-month data sample. The chi- 
square value is 288 kt, which is lowest. The K-S test value is 
also lowest. However, the tally of individual cases ranks this 
method below the others. 

Method B: Up ■ u + F o u , where o u - f(Ri, S, u*). 

The rms error of 3.15 kt is greatest. However, superiority over 
Method A is shown in Table 4-2. 

Method C: Up ■ Gu 

This is a control procedure in the sense that stability is largely 
ignored and a gust factor G(«1.62) is simply applied to u to obtain Up. 
The rms error of 2.81 kt Is smallest and this method ranks highest in 
the tally of individual cases and in Figure 4-7c for 8 u ^ 12 kt. 


The control procedure using a gust factor yielded an rms error of less 
then 3 kt and it was not surpassed by either of two alternative schemes using 
the concepts of stability, the Monin-Obukhov mixing length, and the Richardson 
nunber. 
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